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1 Introduction 



Nonrelativistic quantum mechanics is naturally related to the second order differential equa- 
tions. Exactly solvable quantum mechanics models play a very fundamental role in demon- 
strating the basic concepts and methods of quantum theory. In the textbook fashion, it is 
often the classical results of these differential equations, i.e. the Shrodinger equations, give 
complete solutions to the corresponding quantum mechanics problem. However, in this paper 
we will present two examples where things go in the opposite direction: we can first solve the 
quantum mechanics problems using a combinatorial method, without direct referring to the 
associated differential equations, and then we can say something new about these equations 
based on the combinatorial results. The two ordinary differential equations are the Mathieu 
equation and Lame equation. They are respectively associated to two mathematicians who 
studied them in the 19th century [1], [2]. They arose in classical mechanics problems then, 
and now also arise in simple quantum models. The two equations and their solutions are 
related to each other: the Mathieu equation is a particular limit of the Lame equation. 

The method we use to study the two equations are basically combinatorial, namely we 
will relate the equations to some Young diagrams, and give counting data to each Young 
diagram. With these data we can obtain the eigenvalues following a well defined algorithm. 
This algorithm is rooted in quantum field theory, some beautiful results about four dimen- 
sional M = 2 supersymmetric Yang-Mills(SYM) gauge theory that are developed since 1990s. 
More specifically, we are referring to: 

(1) . Seiberg-Witten theory oiU = 2 SYM[T9l [20], and instanton counting[2T], 

(2) . The correspondence between Af = 2 SYM and classical [30l [311 [32l [33l [341 [35] / quantum [36l 
[37] integrable models. 

This particular class of quantum gauge theory is solved by utilizing its symmetry, duality 
and integrability properties. The intriguing relation between the quantum gauge theory and 
integrable system indicates the two theories share some very basic physical properties despite 
they look very different. In this paper, using this relation, we obtain information about the 
mechanical system from the dual gauge theory. 

The instanton counting deals with Af = 2 quantum gauge theories of general type, with 
various gauge groups and matter multiplets. The gauge theories relevant to the Mathieu 
and Lame equations are: 

(i) The Mathieu equation is related to A/" = 2 pure gauge theory with SU(2) gauge group. 

(ii) The Lame equation is related to A/" = 2 theory with SU(2) gauge group coupled to an 
adjoint matter multiplet. This is the mass deformed maximally supersymmetric Yang-Mills 
theory(A/' = 4 SYM), denoted as Af = 2* theory. 
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More general gauge theories are related to some higher order differential equations, or their 
algebraic form-the Baxter equation. 

The method used to solve the quantum gauge theory /integrable model is in vogue in 
quantum field theory, although is as yet not fully understood. It is hard to understand 
this method solely in the context of quantum mechanics/differential equation, because the 
quantum mechanics Shrodinger equation does not give us any hint about the hidden sym- 
metry and duality properties which are crucial for the solvability of gauge theory. Further 
clarification of the meaning of the combinatorial nature of these differential equations, and 
possible extension to other equations, will be an interesting problem. 

The eigenvalue expansions for the Mathieu equation have been well studied in the liter- 
atures, our combinatorial approach gives exactly the same results and provides some new 
perspective on this equation. For the Lame equation, its expansions are more complicated, 
some analytical results are only obtained recently. We will give an asymptotic expansion for 
the eigenvalue, through the combinatorial method, which coincide with results obtained by 
different method. Some classical literatures containing results of the Mathieu and the Lame 
equations are|3lllEll6l[71IHl[9l[T0llII], several useful original papers are [El [El H 15]. 

2 Differential equations and gauge theory 

The equations we discuss are the 2-particle Schrodinger equations of two kinds of famous 
integrable models: the periodic Toda chain and the elliptic Calogero-Moser model. The 
relation between these classical integrable models and supersymmetric gauge theory is a 
well known result [30| l31] [32] l33t l34t l35] . The quantum gauge theory method can be used 
to obtain some results of the two differential equations, this is based on the following fact: 
According to a recent study due to Nekrasov and ShatashvilifS^ , the vacua of Q deformed 
supersymmetric gauge theories are in one to one correspondence with the spectrum of certain 
quantum integrable systems. The work of Seiberg and Witten[T9l [20] establishes the basic 
notion of four dimensional M = 2 SYM, especially the electric-magnetic duality in the moduli 
space. The instanton counting pro gram [2T] provides a combinatorial approach to the theory. 
Among a wide class of A/" = 2 theories, two theories are simple but typical: the SU(2) pure 
gauge theory [T9] and the M = 2* theory [20]. The two theories are exactly related to the 
Mathieu and Lame equations, respectively, through the gauge theory/quantum integrable 
system correspondence. 

The gauge theories we study carry the simplest nonAbelian gauge group SU(2), the 
corresponding quantum mechanical models are simply two body systems. Remember that 
for an integrable many body system, the building block of its scattering process is the two 
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body scattering data due to the factorization of S-matrix. Therefore, the two body problem 
plays a basic role. After reducing the motion of the center of mass, it reduces the problem 
of a single particle in external potential. 

The gauge /integrable correspondence indicates there is a well defined algorithm to evalu- 
ate the eigenvalues of the Mathieu and Lame equations in a particular asymptotic expansion 
region, up to arbitrary higher order of the expansion parameter, by counting certain Young 
diagrams with weights. With this information we can extend the eigenvalues to other asymp- 
totic expansion regions by utilizing their duality relation. 

2.1 Mathieu equation 

The Mathieu differential equation is: 

-— + {X-2h cos 2^)^ = 0. (1) 
The related modified Mathieu equation is obtained by z — )■ iz: 



{X-2hcosh2z)^ = 0. (2) 



The (modified) Mathieu equation is useful in various mathematics and physics problems. 
That includes the separation of variables for the wave equation in the elliptical coordinates, 
a quantum particle moving in periodic potential. It also appears in problems such as wave 
scattering by D-brane of string theory [17], reheating process in some inflation models p^. 
In the problem we concern, it is the Shrodinger equation of the two body Toda system. 

The potential cos 2z is periodic along the real axes, therefore the Mathieu equation is a 
Floquet(or Floquet-Bloch) differential equation. A Floquet function with period n should 
shift a phase as 

^(z + tt) = e^"^^(z), (3) 

the quantity u = z/(A, h) is the Floquet index. For periodic integrable models, the Floquet 
indexes of wave function are quasimomenta of quasiparticles. Our focus in this paper is the 
eigenvalue A, as a function of u, h. 

The eigenvalue can not be written in elementary functions, only asymptotic expansions 
are obtained when a small expansion parameter is available. And for the Mathieu equation 
there are three asymptotic expansion regions for A, they are located at: 

A ~ z/^ + /^-corrections, h <^ 1, ^ 1 

A ~ ±2h + (^)-corrections, /i > 1. (4) 



4 



These are well known results about the Mathieu equation, they can be found in handbooks 
such as[8l El [TT]. In its relation to gauge theory, each region corresponds to a weak coupling 
description of the gauge theory, they are related by electric-magnetic duality. 



2.2 Lame equation 

The Lame differential equation can be written in several forms. The Jacobi form is: 

■— -[A + n{n~l)kVx]<l> = 0, (5) 

where snx = sn(x|A;^) is the Jacobi sn-function. The n(n — 1) is in accordance with usual 
literature, but we will not discuss whether n is an integer or not which is crucial for the 
classification of the solution. The Weierstrass form is: 



[5 + n(n-l)p(2;)]$ = 0, (6) 

where pi^z) is the Weierstrass elliptic function. The Lame equation is obtained from separa- 
tion of variables for the Laplace equation in the ellipsoidal coordinates. For our interest, it 
is also the Shrodinger equation of two body elliptic Calogero-Moser system. 

The two equations are related by a change of variables [7]. If we change the variables as 

/ = 7 

Vei - 62 

with K' = K{k') the complete elliptic integrals of the first kind and k' = y/l — k^ the 
complementary module. Then Eq. ([6]) changes to Eq. ([5]). We have the relation 

B = {ei-e2)A-e2n{n-l). (8) 

with 



63 - 62 
61-62' 



(9) 



the elliptic module. Therefore 



B 1 + P , , 

A= n(n-l). 10 

61-62 3 

The two equations are equivalent provided the change of variables is smooth. However, the 
change is actually singular in the limit — )■ 0, as in this limit K' 00. Therefore, the two 
equations may reduce to different equations in limits involving /c -t- 0. 

The elliptic functions snx and p{z) are double periodic. The periods of snx are and 
2iK' ] the periods of p{z) are 2a;i and 2uj2. The periods are related by 

K iK' 

OJl = , UJ2 = , (11) 

V 61 - 62 V61 - 62 
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However, note 

sn(x + 2K) = -sn(x), sn(x + iK') = t^—- (12) 

A;sn(x) 

therefore, the periods of the potential sn^x are 2K and 2iK' . 

Thus, the Lame equation also falls into the equations of Floquet type. However, its 
Floquet analysis is not very complete. For example, when the coordinate x shifts a period 
2K, or coordinate z shifts a period 2ui, the phase shift of the function $ is 

$(x + 2K) = e^2^^<l>(x), (13) 



or 

<l>(z + 2a;i) = e^2^^^$(z), (14) 

Where /x and u are the Floquet indexes for equation (jSD and The phase shifts should be 
the same, then we have 

z/ K , , , 

- = — = vei - 62- (15) 

We use different letters [i and v to denote the index for the Lame equation, depending its 
appearance in A or B. In the limit k 0, fi and u coincide. This limit is involved when we 
reduce the Lame equation to the Mathieu equation, and indeed their eigenvalue expansions 
share some common features, see formulae and ( |571) . ( |37j) and ( |82l) . 

The main goal of this paper is to derive the eigenvalue of Lame equation, A (or B), as 
a function of /i (or z/) and n, /c. The Lame eigenvalue also has three asymptotic expansion 
regions, they are located at: 

— {A + — ) ~ /i^ + fc'^-corrections, = n{n — l)fc^ ^ 1, 

— (Ah ) ~ =p 1 corrections, = n(n — 1)/:;^ ^ 1. (16) 

2 2 k 

The second expansion has been given by E. L. Ince, and later by H. J. W. Miiller, see the 
book[TT]. The first expansion was obtained recently by E. Langmann as a special case in his 
algorithim|15]. 

Often, quantities z,X,h and x,A,n{ and z,B,n) are taken to be real, however, in this 
paper we extend them to the complex domain. The reason to do so is that the M = 2 gauge 
theories are related to algebraic integrable systems whose spectral curves are defined over the 
field of complex numbers. Then the abelian variety of the curve, which governs the linearized 
motion of corresponding mechanical system, is the Jacobian variety which carries a complex 
structure, i.e., a complex tori. When we compare the magnitude of these parameters with 
some numbers, we mean taking their absolute values or restoring their real values. 
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2.3 QFT method of deriving the eigenvalue 

If there is a small expansion parameter, then the WKB method is a useful way to derive the 
eigenvalue as an asymptotic expansion in the region where the expansion is valid. Although 
the WKB method is a well studied technique in differential equation theory, we note that 
in the literatures about the Mathieu and Lame equations the WKB analysis is not used to 
derive the eigenvalue. Recently, the higher order WKB analysis was carried out in [38] due 
to the interest about the M = 2/quantum integrable correspondence. The new ingredient is 
that the WKB analysis is carried out in the complex plane, the phase is a contour integral 
along the closed path, therefore a trick of exchanging order of integration and differentiation 
can be applied to obtain higher order expansions. In its relation to gauge theory, the integral 
contours are homology cycles of the Seiberg-Witten curve. Following this development we 
derive all the asymptotic expansions of the Mathieu eigenvalue in|40[ Hi]. However, this 
method needs some higher order differential operators that become very hard to determine 
when the expansion order grows, even though their general structure is conjectured and 
tested up to the 4th order. 

The standard way to solve an integrable system is to derive the eigenvalue of Hamiltoni- 
ans as functions of quasimomenta, and the Bethe equations that constraint quasimomenta. 
The gauge theory /integrable system correspondence, proposed in[36l [37] by Nekrasov and 
Shatashvili, provides a new way to the problem. According to their work, there is a pre- 
cise correspondence between certain quantum integrable models and vacuum sector of gauge 
thoery in the Coulomb branch. The complex scalar condensations {i = 1,2 ■ ■ ■ rankG) in 
the gauge theory are identified with the quasimomenta of integrable model, and the vacuum 
expectation value of gauge invariant operators Tr^'^ are identified with eigenvalues of H^. 
The prepotential of gauge theory is identified with Yang's potential that gives the Bethe 
equation. For G = SU{2) theory, the only nontrivial Tr^'^ is for k = 2, the corresponding 
Hamiltonian H2 leads to the Schrodinger equations. Therefore, for our purpose we need the 
following relation in the gauge theory context, 

< Tr0^ >= u{ai, qin, ei, 62). (17) 

where qin is the instanton expansion parameter, related to some scale of the integrable model, 
ei,e2 are the equivalent rotation parameters. A further limit 62 — )■ is needed, and then ei 
plays the role of Plank constant. 

We noted[l2] that the Matone's relation of A/" = 2 SYM gives a very direct way to 
derive u{ai, qin, ei, £2) from gauge theory partition function. It states a relation between the 
prepotential J-" and m[281 [29] : 

2u = q^n^J^- (18) 
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According to the instanton counting program, the prepotential of = 2 theory can be 
obtained through a combinatorial method. We will show in this paper how to utilize this 
relation to derive the eigenvalues. To fully understand the combinatorial approach, the reader 
should understand some related quantum field theory background, especially the instanton 
counting of A/" = 2 gauge theory [2T| [22], based on techniques developed in [231 1211 125] . 

So first we need the partition function of A/" = 2 gauge theory in the Q background. The 
most nontrivial part is the instanton contribution. Our problem is concerned with gauge 
theories with SU(2) gauge group, the basic combinatorial objects are the Young diagram 
doublet (Yi, Y2). The total number of boxes of the Young diagram doublet is the instanton 
number, |Yi| + 1^21 = k E Z+. For each Young diagram Y^, a = 1, 2, denote the coordinates 
of a box s as {ia,ja), and define its leg-length and arm-length in the diagram Ya,/3 = 1,2, 
as 

hpis) = Pp^i^ - j„, Vp{s) = U,3,jc - ia, (19) 

where i^/j^j^ is the number of boxes in the ia-th row of Yg, while z/^j^ is the number of boxes 
in the ja-th column of Yg. Indexes a and /3 may be not the same, when this happens, some 
boxes s{ia,ja) are outside of the diagram Yg, /^/^(s) and can be negative numbers. See 
Figdll) 
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Figure 1: A Young diagram Y^, with s inside it and s' outside. 

The partition function of the four dimensional A" = 2 SYM theory consists the pertur- 
bative (classical and 1-loop) part and the instanton part, 

00 

Z = ZP^'' ■ Z*""* = ZP""* ■ ^kQin, (20) 

A;=0 

where Qin is the instanton expansion parameter. 

The instanton counting program gives an elegant way to evaluate the instanton partition 
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function. For example, for the SU(2) pure SYM, we have Zq = 1 and for k > l.[26| 



1 



E n 



n E2i(s')(e+-i?2i(s'))i5;22(s')(e+-i?22(s'))' ^^^^ 

s &2 



with 



Eap{s) = - hp{s)ei + {va{s) + l)e2. (22) 

defined by the Young diagram doublets. Here e+ = ei + e2, a^p = aa — ap- For SU(2), we 
have fli = —02 = a. Only some basic doublets are needed, the contribution of (Yii^Ya) is 
obtained from the contribution of {Ya, Yg) by change ^ a^, the contribution of transpose 
doublet (F^, Fj) is obtained from the contribution of (Fq, Yg) by change ei o e2- 
The prepotential J-" of the gauge theory is defined by 

Z = exp(— ) = exp( — (J-^-* + ^ (23) 

The instanton contributions to the prepotential J^k,k = 1,2,--- can be obtained through 
Zk, the first few instanton contributions are 



J^l = 


-eie2^i, 






^2 = 


-£162(^2 


2 




^3 = 


-6162(^3 


— Z1Z2 




^4 = 


-6162(^4 


— Z1Z3 




•^5 = 


-eie2(2'5 


— Z1Z4 


— Z2Z3 4 



^1^2 - ^^1) 



^^f), (24) 



The last ingredient we need is the Matone's relation (1181) . As u also receives instanton 
corrections, u = J2k=o'^k(lin: this leads to 

2no = — ^J-P'^^ 2uk = kJ^k, k^l. (25) 

O In Qin 

As a consequence, if we choose a "good" expansion parameter, denoted as q, then the k-th 
coefficients of g-expansion of the eigenvalue is proportional to the fc-th instanton contribution 
of the SU(2) gauge theory prepotential. The criterion for "good" will be clear when we deal 
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with the Lame equation. This relation is a quantum field theory result, it is fairly novel from 
the point view of differential equation. 

According to the proposal of|36], the quantization of corresponding integrable system 
is achieved when €2 — ?■ 0. Therefore, after we obtain the prepotential, we take the limit 

£17^0,62^0. 

3 Review results of Mathieu equation 

As the Mathieu and Lame equations are closely related to each other, the procedure of deriv- 
ing their eigenvalues is largely parallel. The results in this section have already been given in 
earlier works, see[26l [3H1 SOl HH Il2]. Here we give the results of the Mathieu eigenvalue as an 
educational example to illustrate the basic method and logic of our combinatorial approach. 
After understanding this, it is straightforward to apply this method to the Lame eigenvalue 
problem which is technically more involved. 

3.1 Combinatorial evaluating in the electric region 

First we state a relation between singularities in the moduli space of the gauge theory and 
asymptotic expansion regions of the Mathieu eigenvalue. The moduli space of the M = 2 
SU(2) pure gauge theory is parameterized by the complex scalar field vacuum expectation 
value u = |Tr0^ which breaks SU(2) symmetry to abelian U(l). The physical theory is sin- 
gular at M = ±A^, 00 in the sense that new massless particles appear. These singularities are 
labeled by the U(l) charges of the massless particles, denoted as "electric", "magnetic" and 
"dyonic" , respectively. Near each singularity the gauge theory has an unique weak coupling 
description, therefore perturbative expansion is valid. Our calculation in[10l HI] shows that 
there is a one to one correspondence between the perturbative expansion of gauge theory 
near a singularity in the moduli space and the asymptotic expansion for the corresponding 
Mathieu eigenvalue. For the eigenvalue, in each expansion region the intervals between ad- 
jacent energy levels are very small compared to the eigenvalues themselves, therefore the 
eigenstates are dense there. See Fig([2]) and FigQ. 

The combinatorial method is used to derive eigenvalue at the electric point, at m ~ 00, 
where electric coupling is weak: ^ <^ 1. More explicitly, the parameters of gauge theory and 
Mathieu eigenvalue are identified as 



The exact expression of Z^*^^* is given in[22]. In the limit €2 — )■ 0, we get the perturbative 



A 



8u 



-2 ' 




h 




(26) 
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V(x)=a2cos2x 



electric oo 



netic A/- 




dyonic 



Figure 2: Singularities in the w-plane. Figure 3: Spectrum and duality of cosine- 
potential. 



part of the prepotential[36j: 



.AT N 



i=l 



«J=1 



with Ct7e(x) satisfying 



-^^,(x) = inr(i + -; 

ax e 



(27) 



(2^ 



For our problem N = 2 and = 2a. 

As an exercise and experiment, we work out the h^^ order of A for h 1, i.e. the 

five instanton contribution. The instanton counting formula is ( |2T]) . For the 1-instanton 
contribution, we need only one Young diagram doublet: 

(D 0) 

where is the empty diagram. Then we get 



eie2(a?2 - (^i + £2)^)' 
For the 2-instanton contribution, we need two Young diagram doublets: 



(29) 



and we get 



2a?2 



- 17eie2 - 8el 



t\t\{a\2 



(ei + e2)2)(a?2 - (2ei + e2Y){al2 - (ei + 2e, 



(30) 



Similarly, we need three Young diagram doublets to derive Z^: 

. 0) ( rm . 0) (cn 



□) 
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^3 

num 
den 



numerator 
denominator 

-2(20^2 - 1<oa\^e\ + - ^la\^exe2 + 363e?e2 - 260^2^2 
-594e?e^ 

3e?e^(a?2 - (ei + e2)')(a?2 - (ei + le^f^al^ - (ei + 362)^) 
(a22-(2ei + e2)')(«?2-(3ei + e2)'), 



+594e?e^ + 363eie^ + 72e^), 



(31) 



For Z4 there are seven Young diagram doublets: 



, 0) 







□ 



□) 







Z4 = 
num = 



numerator 



J.3 

1^2 
^2 



denominator 

4a^2 - 132a^2e? + 14400^2^^ - 6208a?2e? + 9216e? 
-212a^2eie2 + 56440^26^2 - 44336a?2e^e2 + 100608eje2 - 1320^2^2 
+8651a^2e?e2 - 1241390^2^^2 + 440688e?e^ + 5644a^2eie2 - 171845a?2e?( 
+1009131e^e^ + 14400^2^ - 1241390^26^6^ + 13199946^6^ - 443360^2616^ 
+100913l6?6^ - 62080^2^2 + 4406886^6^ + 1006086i6^ + 92166^, 
den = 66^6^(0^2 -(ei + e2)')(a?2-(ei + 262)')(a?2-(ei + 362)2) 
{a\^ - (61 + 462)')(a?2 - (261 + e2)')(a?2 " (3ei + 62)') 

(a?2 - (461 + e2)')(aL - (26i + 262)'), (32) 

With Zfc in the hand, it is easy to derive J-'k using f l2^ . The complete form of the first 
three ^-"1,^-2,^-3 are presented in[2S]- As observed in|12], according to the Matone's relation 
( 125|) , and taking in to account (126|) . the Mathieu eigenvalue is 



00 ^4fc-2 

A = z/2 + ^^-p4fcJ'fc(ai2,6i,62 = 0)/i 



16* 

k=l 

00 , 
2 _ , h 



+ ^4A;J-,(z/,l,0)(-)2^ (33) 



k=l 

,2 



The leading term z/ comes from the perturbative contribution. We have used the fact that 
the numerator and denominator of ^-"^(012, 61) are homogeneous polynomials of degree 2k + 2 
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and 6k, respectively. Up to A; = 4, we get 

2 1 j,2 5Z/2 + 7 4 9z/4 + 58z/2 + 29 
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2(z/2-l) 32(z/2-l)3(iy2_4) 64(i/2 - 1)5(^/2 _ 4) 

1469i/^° + 9144Z/S - 140354Z/6 + 64228z/^ + 827565z/2 + 274748 , o 



h^ + --- (34) 



■2-9) 



8192(z/2 - l)7(z/2 -4)3(zy2 _9)(zy2 _ ^ ^ 

for /i < 1, 1/ > 1 and/or /i > 1, ^ < 1. 

3.2 Extension to other regions 

A complete solution of the eigenvalue problem requires extending the asymptotic expansion 
from M ~ oo to other two regions near u = ±A2. This is achieved by utilizing the duality 
relation between the three regions, a crucial result of Seiberg-Witten theory. Physical quan- 
tities at u ~ oo have their dual counterparts at m = ±A2. For undeformed gauge theory, 
ei = €2 = 0, the scalar condensation a and its dual are given by integrals along the 
homology circles of the Seiberg-Witten curve which is a torus for SU(2) gauge theory. As 
functions of u, they can be expanded in the three regions. The leading order WKB expansion 
of the Mathieu equation gives solution of the undeformed gauge theory [381 1391 HO]. 

Theory with ei 7^ 0, e2 = is the "quantized" version of theory with ei = €2 = 0. As first 
studied in[38], the WKB method, widely used in quantum mechanics problems, provides an 
efficient method to obtain higher order quantum corrections. Actually, the work of [38] shows 
that the ei-corrected scalar condensation and its dual can be generated from the classical 
one by action of certain differential operators with respect to u and A. 



where \sw is the Seiberg-Witten form, in the context of gauge theory /integrable system, 
the integral / Xsw coincides with the classical action J p{z)dz of a particle. D{u,du,ei) can 
be expanded as D{u,du,ei) = 1 -|- J2'^=i^"'^n{u,du)- The differential operators Dn{u,du) 
are polynomials of u and du, the first few of them can be found in [381 SOI SI]- Reverse the 
function a(u,gj„,ei) we can get u(a,gj„,ei) which gives the eigenvalue through parameters 
identification ( 126|) . 

The WKB quantization method is independent of the instanton counting approach, in 
principle it serves as a self-consistent way to solve the eigenvalue problem. Nevertheless, the 
WKB method gives a nontrivial check for the combinatorial method. It turns out that at 
u ~ 00, series expansion of a{u, A, ei) from (!35l) exactly gives eigenvalue ( 1341) . 

The dual an is given by integral along the dual homology circle. 




(35) 




(36) 
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D{u, du, ei) remains the same, and a reverse of the function aD^u, qin, ei) to u = ei) 
is needed. Near u = where the magnetic couphng is small, ^ ^ 1, we found that the 
series expansion of a^i gives another well known eigenvalue expansion: 

X = 2h- AuVh H \ ^ H — 

528z/^ - 1640Z/3 + 405z/ , , 

+ , + • • • 37 

for /i 1, ^ ^ 1, with identification 

, 8m 2iaD , 4A2 4^^" 

A = -J, z/ = ^— , /i = — = — (38) 
ei ei 

The details are given in[lT], see there for result up to the order h~2. Moreover, the third 
expansion region near u = — A^, the dyonic region, is the mirror of region near m = A^. The 
eigenvalue there is obtained through u ^ iu, h ^ — /j,[40]: 

4z/2 + 1 4i/3 + Su 80z/4 + 1361^2 ^ g 



A = -2h + AuVh 



23 2^Vh 212/i 

528i/^ + 1640Z/3 + 405z/ 



216/if 



+ ■■■ (39) 



We should stress that we can combine the combinatorial and WKB method to make a 
calculable algorithm to get all eigenvalues to higher order. As we have said, the higher order 
differential operators used in[38lllT] become very hard to determine. With the combinatorial 
results of higher order electric expansion, according to the Claim 1 of^Tj, we can easily 
extend the WKB differential operators to even higher order. It becomes quite straightforward 
if some simple computer codes are used properly. For example, only a bit more labor is needed 
to get the next order results, continuing |41j: 

4471Z/12 + 69361i/i° - 1039598//*^ - 28444301/^ + 13541915z/^ + 20651309z/2 + 4453452 

16384(z/2 - l)9(z/2 - 4)3(z/2 - 9)(z/2 - 16)(z/2 - 25) ' 

(40) 

for /i <^ 1, ^ 1 and/or /i 1, ^ <^ 1. Now we can determine -D„(m, du) up to ri = 10, and 
then continue to determine higher order magnetic expansion. With this strategy, we get 

^^(21635328i/^° - 388041984//^ + 1537405408i/^ - 1708613456z/^ + 449869257//^ 

- 11917692) + ^^(4295427072//" - 97098594560//^ + 523957083264//^ 

- 888479287968//^ + 440967453876//^ - 41540033277//). (41) 

for /i > 1, ^ > 1. 
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4 Combinatorial approach to Lame eigenvalue 



The Lame equation is closely related to the Mathieu equation. In the limit n — )■ oo, — )■ 
while keep n{n — l)k'^ = fixed, the Jacobi sn-function becomes the trigonometric 
function sin, and Eq.(l5]) becomes equivalent to Eq.([l]). This fact indicates each asymptotic 
expansion of the Mathieu eigenvalue presented above should be the limit of a Lame eigenvalue 
expansion. That means the Lame eigenvalue has at least three asymptotic expansion regions, 
located approximately at three point as in ( |T6|) . In this section, based on the quantum field 
theory, we argue that there are only three asymptotic expansion regions. We also derive the 
eigenvalue at the electric region using the combinatorial method, and write it in a proper 
form. 

4.1 Identify parameters 

First, we have to make a identification between the module parameters in the Lame equation 
and in the M = 2* gauge theory. The parameter k appearing in Eq.(l5]) is the elliptic module 
associated to the periods AK{k) + 2iK' (k) of the Jacobi elliptic function snx. The Jacobi 
theta constants 6i are often expanded in another module parameter, the nome q. The two 
parameters are related to each other, their precise relation can be found in Appendix |8l 
Especially, we notice 

1 3 5 9 (a)'^ 

e = 16g5 - I28g + 704gi - 3072g2 + 11488g5 - 38400g^ + ■ ■ ■ = -^r^- (42) 

Our claim is that q is exactly the instanton expansion parameter of the M = 2* theory, and 
is the "good" expansion parameter for the eigenvalue B in the electric region. 

Another thing we have to identify is the relation between the parameter n in the Lame 
equation and the adjoint mass of the gauge theory. The equivalent parameter m appears in 
the instanton partition function of A/" = 2* theory as|27]. 

As clarified in|13], the parameter m differs from the physical mass m* by 61,62 shift: m = 
m* + (61 + 62)72. We identify n = e^^m. 

Lastly, we have to identify the relation between the eigenvalue A, or B, and the scalar 
condensation u. The Seiberg-Witten curve of the M = 2* theory is 

y"^ = {x — eiu e\w?){x — el'm'^){x — e^u elni^), (44) 
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withlfl 

u =< Tr02 >=u + m'^i^ + Ciq + Cag^ + ■■■)• (45) 



We have rescaled the mass in [20] by m — )■ ^ in order to match the mass parameter that 
appears in the instanton counting formula (H3|) . Again this is an elhptic curve, with two 
conjugate homology circles a and /3, and the Seiberg-Witten differential one form, u and u 
differ by terms caused by mass deformation. It turns out that u is directly related to the 
eigenvalue B. 

In summary, we identify 

n — —, JJ — ^ — ' 1 — lin- I4bj 

ei ef ei 




Figure 4: Singularities in the w-plane for N = 2* theory. 
The moduli space has three singularities in the u plane, as in Fig.(j4]), they are located at 

it = -^eim^ = m2(-^ + 2g + 2g2 H 

o 12 

u = ie2m2 = m2(-^ - g5 _ g _ 4gf _ g2 _j 

o 24 

_ 1 1 1 3 

u = —e^m^ = m?{ V q"^ — q + 4g2 — . . (^47) 

8 24 

Note that the stationary points of the potential p{z) are exactly at p{z*) = 61,62,63 where 
dzp{z) = 0,d'^p{z) > 0. In the u plane, the singularities are located at 

u = m2(^-4g2-l2g^ + 12g^ + ■■■), (48) 
8 

u = =Fm2(^-3g + 4gi -Tg^H ). (49) 



^Thc relation proposed in^20jis u =< Tr^^ >= u + j^eivn? . But as noted in[44], the one instanton 
calculation for the coefficient ci indicates this relation is not correct. The correct relation is pointed out in 
[45] as u = u — ^ + ^E2{t), with E2{t) the second Eisenstein series. Some scaling of parameters in that 
paper are needed to match ours. 
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As in the case for the pure SU(2) SYM(Mathieu equation), there are only three singularities 
in the moduli space, therefore, there are only three asymptotic expansion regions for the 
eigenvalue. While the Mathieu equation has two free parameters h and z/, the Lame equation 
has three free parameters, n, k and z/, it is not a prior fact the number of expansion regions 
remains to be three. However, the physical argument indicates there is no other expansion 
region. 

4.2 Combinatorial evaluating in the electric region 

We derive the eigenvalue near the electric point. Let us look at the perturbative part, in the 
limit e2 = 0, the prepotential is given by [36] 

N N 

-ppert _ ^^^^~^^2 _|_ 2Tiie^ (tu^^ (fljj ) — w^^{aij — m — ei)). (50) 

i=l *ii = l 

Here we have N = 2, and = exp(27rir). Only the first term a^lngj„ contributes to the 
eigenvalue. 

Then we come to the instanton part of the gauge theory, or the q correction part of 
the Lame equation. Following the line of instanton counting, we derive the prepotential 
^fc(ai2, €i,e2,m),k = 1, 2, 3, ■ ■ ■ , and set €2 = 0, the eigenvalue B can be obtained via the 
Matone's relation. The whole procedure is the same as that for the Mathieu eigenvalue, but 
with new counting formula (H3l) and identification (H6!) . 

Without presenting the results of Zk, k = 1,2,3 ■■■ , which are fairly lengthy, we only give 
the final results of the prepotential. The 1-instanton contribution is 

2[m'^ — 2m^ei — m?{a\2 — 2e\) + m(a^2ei — e\)] 

J~ iQin 2 2 

— Ci 

_ 2m(m - ei)[m(m - ei) - aj^ + ej] 

2 / 



a 
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Writing in terms of n, v, q, it is 

_ r,2n(n — l)q r / n / •> ^ / n 

J'lq^n = el \ ^ X [n{n - 1) - {v^ - 1) , (52) 

Similarly , the 2-instanton contribution gives 



-I2n\n-l)\u''-lf 
+Qn{n - l)(z/2 - lf{u'^ - 2) 

-3(z/2 - l)3(zy2 _ 4)] (53) 
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The 3-instanton contribution gives 

- -■ 3(.'^ - IH.-"- it' - 9) ' P"°<" - + ''"^ + 

-8n^(n-l)^(z/2-l)2(7z/2 + 17) 

+4n^(r2 - l)3(z/2 - l)2(15z/^ - 37u^ - 2) 
-24n2(n-l)2(z/2-l)^(z/2-4) 
+3n{n - l)(z/2 - l)4(z/2 - - 4) 

-(z/2-l)5(zy2_4)(^2_g^j ^54) 

A notable fact is that in the final results of instanton counting the parameter n always 
appears in the form n{n — 1), consistent with its appearance in the Lame equation. This 
is already visible from the numerator of counting formula where terms involving m 
combine to m(m — ei). The full prepotential is 



oo 



^ = J^P'"-' + ^ J^kqt = T^-^' + J-ig + J-2g2 ^ T^i ^ ■ ■ ■ (55) 



fc=i 



According to the Matone's relation, the scalar condensation u is given by 1u = Q^^, 
and according to[45j, u and u are related by 

, m(m-ei) m(m - 
w = M + E2 (r) . (56) 

where we have turned on the ei (quantum) correction to the "classical" relation, and £'2(1") 
is the second Eisenstein series. Now we are in the position to get the eigenvalue B: 

^ = = -1!%-^- ^1-^(1 -^^'^M) 



2 n{n-l) 8n{n-l)q / 2 im 

/2 - ^ ^ (1 - 2^2(9)) \ -^[n{n - 1) - (z/2 - 1) 

6 I/^ — 1 



X [^=^(72 - l)3(5z/2 + 7) - 12n\n - - 1)^ 



(z/2 - 1)3(1,2 _ 4) 

+6n(n - l)(z/2 - l)2(z/2 _ 2) - 3{u^ - - 4)] 

-8n^(n - l)^(z/2 - l)2(7z/' + 17) 
+An^{n - l)^(z/^ - l)^(15z/^ - 37//^ - 2) 

-24n2(n - l)2(z/2 - l)4(;y2 _ 4) 

+ 3n(n - l)(z/2 - l)4(^y2 _ 3)(^2 _ 4) 

-(i/2-l)5(zy2_4)(^2_g)]^... ^57) 
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We have used the perturvative result qdqj^^'^^^ = efu'^/A. In the electric region, we have the 
leading order behavior 8u ~ m^, that is z/^ ~ n{n — 1). Therefore, B ~ —^n{n — 1) + 0{q). 
The expansion of B is valid when 

g < 1, n ~ z/ > 1, ng3 < 1, (58) 

and it is also valid for 

g < 1, n ~ z/ > 1, 1 < ng3 < z/. (59) 

An equivalent expansion was also obtained by E. Langmann in [15]. His method is differrent 
from ours, we have shown the equivalence of his eigenvalue expansion S and our expansion 
B, see Appendix [TOl 

In order to compare to the WKB results, we do some further works to derive the expansion 
of A from B. B can be further expanded as 

B = -z/^ + -72(72 -l)(l-24g-72g2-96g^ ) 

3 

+ 6g' + 12g=^ + ■ ■ ■ ) 



2 1 

-z/^ + -n{ 


[n - 


o 




8n^{n — 


I? 


z/2 




8n^(n — 


I? 


z/4 




8n'^{n — 


I? 


z/6 



[q + (18 - 12n{n - l))q'^ + (84 - 96n{n - 1))^^ + ■ ■ ■ ] 

■[q + (66 - 60n{n - 1) + 5n'^{n - 1)^)^^ 
+ (732 - 960n(n - 1) + 2A0n^{n - lf)q^ + ■ ■ ■] 

+ ■■■ (60) 

Then it is easy to get the expansion for eigenvalue A, taking into account the relation between 
A and B ffTOj) . and the relation between u and fi f|T5|) . As we have worked out terms of order 
up to q^, taking into account relation we keep terms up to the order fi~^ and k^'^. We 
finally get 

A = -^^--n(n-l)e-n(n-l)(—k' + —k'^ + ^k'' + ^k'^ + ^^k'^ + ---) 
^ 2 ^ ^ ^ ^^6 32 2048 4096 65536 ^ 



/i2 ^32 4096 4096 131072 
n\n - 1)2 _ _ ^ 7 + 6n(n-l) ^ ,^s _ ^ 7 + 6n(n - 1) ^ ^^^p 



g4 V 4ogg ^ V 8192 



.73 + 60n{n- 1) 



131072 

^2(^-1)2 1 . 1 ,74 - 6072(72- 1) +5722(n- 1)\,„ 

^ ^ f -A:'^ r + ( ^ ^ —)k^ 



fx^ '32 32 ' 8192 

, 29-10222(72-1)2 ^,^,, , 1 

^ 131072 

+ ••• (61) 
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Before doing the WKB check, we first give a simple consistent check using a hmit of 
the Lame equation. The M = 2* theory has a decouphng hmit. When the energy scale is 
very small, or the mass is very large, the adjoint mass decouples and the theory flows to the 
pure gauge theory. In this limit, every singularity in the moduli space of the M = 2* theory 
becomes a singularity of the pure gauge theory. The limit is 

g 0, m cx), m^^g-^-A^, (62) 

or 

^ 0, m ^ oo, m^k'^ -16A^. (63) 
Writing in the parameters of the differential equation, the limit is 

h 

g — )■ 0, n — )• oo, n(n — l)yg — )■ — — , (64) 

or 

A; 0, n oo, n{n - l)k'^ -Ah. (65) 

Accordingly, in this limit the Lame equation reduces to the Mathieu equation, and the Lame 
eigenvalue ( 157|) and ( 161]) reduces to the Mathieu eigenvalue ( !34|l : 

~B^X, -A^\-2h. (66) 



5 Check through the WKB method 

When there is a small expansion parameter, the WKB method can be applied. This is 
well studied technique in some differential equations, as we have said, the new ingredient 
contributed in [38] is doing the WKB analysis in the complex domain. In this section, we do 
the same thing for the Lame equation, as have done for the Mathieu equation in|38l HOI HI]. 
Rewriting the Eq.(l5]) as 

- [cj + sn^x]* = 0, (67) 

where 

2 2 A 

n(n — 1)A;^' n{n — l)k'^ ^ ^ 

Note that ei is contained in e and to. Suppose e <^ 1, i.e., nq^ ^ 1, then e is a valid small 
WKB expansion parameter. We choose Eq.(l5]) to do the WKB analysis because the readily 
available integrals formulae are suitable to deal with the Jacobi elliptic functions. It is legal 
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to choose Eq. ([6]) instead, as in[50l |5T] where the leading order expansion is performed. 
Expanding the function $ as WKB series: 



= expi / dx p[k 



expi I rfx'(^^^ + pi(x') + ep2(x') + e2p3(x') + ■■■), (69) 



then we have 

Po = iV2Vuj + sn^x, pi = ^(In(po))', 

1 Q ' " 

Apo 2 Po Po 2 Po 

P4 = -(— ) - — , ■■■ (70) 
2 Po 2po 

where the prime denotes 

We have to evaluate the contour integral i^^podx, where a is the homology cycle of the 
curve dH]) relevant to electric scalar condensation. Introduce the amplitude, 



if = amx, (71) 

then we have 



snx = sin(y9, dx = — ^ (72) 

\/l — k'^ sin^ Lf 



The integral becomes 



u) Vu) + sn^xdx = 2 / Vuj + sn^xrfx = 2 / W ^ "^^^^ 'f d^p^ (^73) 
/a Jo Jo \ ^- sm 9? 

In the electric region, A ^ h?, therefore cu ^ 1, we have expansion 

I ~y sin^o9 sin'^09 sin^o? Ssin^o? , 

All the integrals we need to do are of the form 



- • 2r) 

2 sm 



^ di^, with A = y 1 — sin^ n = 0, 1, 2, 3, ■ ■ ■ (75) 
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They are given in Chapter 2 of [TO]. After collecting all terms together, we getO 



I I ^, or. ^ 1/2 (2 + A;2)K - 2(1 + P)E 

^ + sn2xrfx = 2irv^ + -^^^ ' - ^^^4 ^ 



(8 + 3fc^ + 4fc^)ir - (8 + 7A:^ + %k^)E 

(48 + 16P + I7k^ + 24A;6)i\: - 8(6 + 5P + 5^4 ^ g^e^^ 

1344fc8 

+ •••• (76) 

where K and are the complete elliptic integrals of the first and second kind, respectively. 
For higher order contour integrals, as usual the odd order integrals do not contribute, 

P2i+id>i = Q, 2 = 0,1,2--- (77) 

Hence, all odd differential operators D2i+i{u}, d^, k) = 0. While for the even order integrals, 
using the trick of [38], we get the operator D2{oj, d^, k) from 

/ i sn^xcn^xdn^x 

a,/3 Ja,l3 0\/ Z [u + SU^x) 2 



+ {1 + e + 3ue)d^ - h^] (f) podx. (7J 

4 Ja.B 



■^Instead, we can also expand 



1 H sm (p H sm (p H — sm ip H sm (p - 



- sin2 ^ 2^8 ^ 16 ^ 128 

The integrals we need to do become 

/ ' Asin^" ipdLp, with A = Jl-P sin^ ip, P = - -. 
Jo ^ ^ 

The formulae are also given in Chapter 2 of [TU]. This gives the same final result (|7S)) . 
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and get D^^u, d^j, k) from 



a, 13 



1 / G{x) 5 (snxcnxdnx)^ 




^ r ,7 (G(x))2 1 G"(x) 




where G(x) = 1 - 2(1 + A;2)sn2x + 3Pgj^4^^ ^^^^^ _ ^^^(x). 

The contour integral along the a circle is actually the phase shift of the function $(x) 
in ([T3D, i.e. 



The remaining work is straightforward but tedious: expand the po integral fl76l) as power 
series of u and fc, use the differential operators from fl75]) and fl7^ to generate the next two 
integrals for P2,P4, we get the series expansion of /i = /i(c<;, A;). Then inverse the series /i(co', fc) 
to get a series w = uj{fi,k), we finally get the series expansion for the eigenvalue A using 
A = n{n — l)k'^u. We have checked it indeed gives the expansion ( 16 ip . up to order and 
k^. In our WKB analysis, the fc-expansion is truncated at order k^ because in ( I76p we have 
only worked out terms up to the order o;"'^/^, with the accuracy of order k^. At this order 
the WKB check matches sixteen coefficients appearing in ( 16T]) . 

6 Extension to other expansion regions 
6.1 Magnetic expansion 

The last task we have to do is extending the asjnnptotic expansion in the electric region to 
the magnetic and dyonic regions, as what we have done for the Mathieu equation in |?T]. 
We stress that the asymptotic expansion in the magnetic region has been worked out in the 




(80) 



Therefore the Floquet index is given by 




(81) 
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literature, using purely mathematical technique. The explicit formula is given in the book 
by H. J. W. Miiller-Kirsten pT| . which cites results from earlier original papers of E. L. Ince 
and H. J. W. Miiller. It is also presented in Chapter 29 of [9]. It is 

A = -22K/i-^(l + fc2)(V-l) 

' [(1 + P)2(4^3_3^) _4P(4^3_5^)] 



^ [(1 + k^Y{528fi^ - 1640/i^ + 405/i) - 24A;2(1 + k'^f{n2fx^ - 360/i^ + 95/i) 



2l3fi;3 

+16A;^(144/i^ - 520/i^ + 173/i)] 

+ ••• (82) 
The expansion is valid for 

nk > 1, nA; > fi. (83) 

Note k r-u qj_ We remind the readers that the Floquet index in[TT] is denoted by q, its 
relation to /i is g = 2i^. Also their eigenvalue is denoted by A, related to ours by A = —A. 
Again we use the letter fi to denote the Floquet index of A. In the limit n — )■ oo. A; — )■ 
while keep k"^ = n{n — l)fc^ — —4h fixed, the eigenvalue fl82|) reduces to the corresponding 
Mathieu eigenvalue, see formula ( 137|) . 

Now, let us look at how to extend the Lame eigenvalue (1611) in the electric region to the 
magnetic region. The tool is the WKB data obtained in the last section, and apply them to 
the /3-contour integral. In the magnetic region located near A ~ 0, where w <^ 1, we can 
expand v^oT+sin^ as: 



2 u 1 1 u"^ 1 5u^ 1 

2 simp 8 sm Lf losm 128 sm 

All the integrals now we need to do are 



1 



Asin^"-^ Asin^^-V 



d^ = '^ . . 2n-l (85) 



with 

1 



sinyjQ = -r, ^ = \l 1 — k'^ sm. (p, ?t, = 0,1,2,3, 



They are also given in Chapter 2 of [10]. The first few orders of the integral '^^^podx give 

r , ^, .A 1 + A;2 3A;^ + 2A;2 + 3 o 5(1 + A;2)(15A;4 - 6A;2 + 15) , 

(7) Vo; + sv?HdK = mi-uj cu^H cu" ^ 

fs ^2 16 128 6144 



^7) 
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Using the differential operators obtained in ( ITHj) and ( !79|) we liave obtained in tlie last section, 



There is a point a bit puzzling to us. When we go through the circle /3, the coordinate x 



index jj,. But we find the correct number to divide by is iir. 

The remaining work is the same as that in the last section: use the differential operators 
in ( 1781) and (179|) to generate the next two integrals about P2,P4, so we get the series expansion 
of fi = fi{u, k). Then inverse the series /i(a;, k) to get a series u = uj{fi, k). Finally we get the 
series expansion for the eigenvalue A using A = n{n — l)k'^uj = n'^u. Or we can substitute 
(!82|) into the series fi = fi{uj, k), and check if the relation holds. We have checked the first 
few terms in (l82l) . up to terms of order k"^. 

From the magnetic expansion A{fi, k) given above, using the relation ( fTOl) and changing 
/i, k to z/, q, we get the magnetic expansion of -B(i^, q). Of course, -B(z/, q) is not economic for 
the magnetic expansion. 

6.2 Dyonic expansion 

The eigenvalues near the dyonic point are mirror to the ones at the magnetic point. This is 
most obvious for B{v,q) ~ u. From ( H71) we know that the magnetic point u = 62^^/8 is 
mapped to the dyonic point u = e^m? by 62 — )■ 63, or — ). — ga. Also from the story of 
Mathieu equationjlO], we know that the index of magnetic region is mapped to the index of 
dyonic region by — ?■ iv. Therefore by the mirror map q^ — t- — p ^ iu we get the dyonic 
expansion Bd from its magnetic expansion Em- 

It is also interesting to look at how the dyonic expansion of A is mapped from its magnetic 
expansion. For this purpose, we need to know how the other parameters are changed by the 
map. It simple to see that under q^ — — g^, we have 62 and 63 interchanged while ei 
unchanged. Furthermore, since 



we can generate integrals p2d>c and p^dx. 

For the Floquet index in the magnetic region, our claim is 




(88) 



shifts 2iK' = 2iK(k'). Therefore the phase (fiapdx should be divided by 2iK' to get the 



fc = 4g3 J]( 



00 



l + q 



,n 



4 



* = n( 



00 



1 - g" 2 4 



(89) 



n=l 



1 + g"-2 



n=l 



1 + g"-2 



The map leads to a simple change for k and k' 





(90) 
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As for the magnetic expansion, A and B satisfy 

Bm = {ei - e2)Am{fi,k) - e2n{n - 1) 



= iei-e2)Ami , " ,k{q))-e2n{n-l), (91) 

The mirror map changes it to 

Bd = (ei - e3)A„(^i^=, i^) - - !)• (92) 

Vei - eg k 

Then substituent this into the dyonic relation B^ = (ei — e2)Ad{fi, k) — e2n{n — 1), we get 

Ad(/i, A;) = — — —Am{—?^^^=,A)-— — —n{n-l) 
ei - 62 vei - es k e\- e2 

ei-eg /ei - 62 .A; 63-62 , - , . 

-A[i^\\ , I—) n(n - 1). (93) 



61-62 V ^1 - 63 k! 61 - 62 
where the function A is the magnetic expansion fl82|) given in the previous subsection. Using 

= k'\ = k\ (94) 

61 - 62 61 - 62 

it simphfies to 

Ad = k''A{^,'^)-kMn-l), (95) 
with the first few terms expanded as 

Ad = -K^ + t2Kfi + ^{l-2e){^ + l) 

Zj rL rv rv rv rv hj 

A WKB analysis can be carried out for this expansion. Set A = —k^ + A, i.e., u = 
then the Lame equation fl^ becomes 

-[u- cn^x]^ = 0, (97) 

The WKB solution for this equation is very similar as in the magnetic case, now with the 
leading order solution po = i\/2\/oj — cn^x, and the integral contour is 7 = — (a + /3), 



{iV2y^ (h podx = -2 - cn^xc/x = -2 J ^d(f. (98) 

J-y Jo Jo V ~ ^ ^ 

In the dyonic region, the eigenvalue is expanded for k ^ 1,k ^ /i, i.e. ^ 1, therefore 
using 

nz 7, — / bJ \ \ Zt" \ 5cD^ 1 N / N 

^UJ - COS^ LP = Z COS (p - — — — — — + ■ ■ ■ , 99 

2cos(p 8 cos-^ ip lo cos'' 128 cos' 
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the integrals we need to do are 



1 /"^° 1 

T-d^ = ^ A.r..2n-iJ V^ n = 0,l,2,3,---. (100) 



A cos^"~^ (p Jq A cos^*^"^ (p 

The first few of them gives 

(bVuj- cn^xfix = -Tc(—u-\ rr^-w H tt^ w H — ^ 777; -uj*+- 

(101) 

For the higher order WKB integrals, as in the previous cases, for odd order integrals 
<^^P2i+idx = 0. For even order integrals, they are generated from (^^p^dx by some differential 
operators. We get D2{co,d^,k) from 

_ , f i sn^xcn^xdn^x 
P2a>c = (J) — -= — —d>c 



(u — cn^xl 2 



8^2 

._L[(_1 + k^ + 2^- Ak^u + 3k^^^)dl 



+{l~2k'^ + ?,k'^i:j)d^-^k'^]j)PQdx. (102) 



and D4{u,d^,k) from 



~w - / 1 r7 (G'(x))^ 1 G"(x) 

^ ~ f^64v^2^12(25-cn2x)V2 12 (^5 - cn^x)^^^^^ 

= 4{7^[(1 - '^')' - 4(1 - 3^ + 2k^)u + 2(2 - 11A;2 + llA;^)^' 
64 45 

+12fc2(l-2A;2)w3 + 9A;4w^]9i 

-^[(1 - 3k^ + 2A;4) - 2(1 - 5P + 5A;^)w - 8k\l - 2^)^^ _ 8A;4^3]^3 
+ i[(10 - 43A;2 + 43A;^) + 66k\l - 2^)^ + ggk^'u^p^ 



+2[k\l - 2^) + 3A;^w]95 - h^}(^^pod>c. 
We find that for the dyonic expansion, the Floquet index is given by 



(103) 



^ = — (j)dxp{x). (104) 

At the moment we do not have a satisfying explanation for the factor fc /vr, but the relation 
is indeed satisfied. We check this up to terms of order k"^ in fl95|) . 

However, some interesting phenomena are observed. Look at the leading order contour 



integrals for the magnetic and dyonic case, formulae ( jHTl) and fllOll) . they are related in a 
simple way: 

k — )■ i-7, uj — — cD, leads to: - (f) a/cj + sn^xrfx — )• — (f) — cv?>cdx. (105) 
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Note that e ~ k ^ ~ ^ also changes as e — )■ —ik'e. Then look at the differential operators, 
formulae (1781179^ and fll02lll03p . they are also related in a simple way: 

k ^ ~ _ 

k^i-j, w — )■ — w, leads to: e^Dn{<jJ,d^,k) ^ e^DJijO^d^^k)^ n = 1,2. 

(106) 

Combining the two facts together, under the mirror map we have 

dHpix) — i- — (J) dxpix), i.e. /x — )■ i^. (107) 

This is consistent with ( !95l) where the index appears as i/i/A;'. 

7 Conclusion 

Combinatorics is a very interesting and fruitful subject of mathematics, its relation to inte- 
grable models is not new|46]. In recent years, some progresses continue to reveal its fascinat- 
ing connection to quantum field theory, integrable model, and string/M theory [22l [36l Wl\ HS] . 
Based on a relation between Af = 2 theory and quantum integrable system[36l EZ], we pro- 
vide a quantum field theory approach to the eigenvalues of Mathieu and Lame equations. 
The approach is combinatorial, based on the instanton counting program^]. 

Our main result is the expansion formula (157|) and (1951) . At the moment we know some 
other independent methods that can derive the same results. The first one is E. Langmann's 
results [I5] that give an expansion equivalent to ( 157|) . we have shown the equivalence in 
Appendix [TOl The second is the WKB analysis we give in section 5 and 6, giving all three 
expansions. Also in the decoupling limit the expansions correctly reduce to expansion of the 
Mathieu eigenvalues. 

We can achieve a relatively simple expansion formula of ( l57l) partly because we have cho- 
sen q as the right expansion parameter. A similar lesson comes form the recently discovered 
AGT relation(Alday-Gaiotto-Tachikawa) |48] . which connects four dimensional M = 2 gauge 
theory to two dimensional Liouville type conformal field theory (CFT). There, expanding the 
CFT conformal block according to the gauge theory partition function has shielded new light 
on nonrational CFTs. On the other hand, studying the Lame equation of the form ([5]) in 
the 6th section reveals interesting duality properties of the N = 2* theory. 

In this paper we do not discuss the combinatorial connection to the eigenfunctions: the 
Mathieu function and Lame function/polynomial. There are examples where a full combi- 
natorial approach to the eigenvalue and eigenstate exists. The Calogero-Sutherland model 
is a famous example, see a very recent discussion in [52] , and relevant references therein. In 
fact, the Calogero-Sutherland model is the trigonometric limit of the elliptic Calogero-Moser 
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model. In another recent paper [19], it is argued that the AGT relation with full surface oper- 
ator, in the semi classical limit e2 — ?■ 0, gives the eigenfunction of the corresponding quantized 
Hamiltonian. See also|51j. This maybe helpful for an explicit combinatorial construction of 
the eigenfunction. 

While the basic tools of the paper, such as Seiberg-Witten theory, instanton calculation, 
have been extensively studied by mathematical physicists and have a solid mathematical 
foundation, the knowledge in this field is still expanding. Some points of our results lack 
rigorous mathematical proof. This includes the general form of -D„(u, du) conjectured in[3T], 
the puzzling factors {iir)^^ and k'(Tc)~^, the convergence of the eigenvalue expansions, the 
precise meaning of the potential J-" for the equations, the relation to Langmann's method. 
We hope to clarify some of them in the future. Therefore, we do not present a mathematically 
rigorous treatment of the Mathieu and Lame equation in this single paper, but just to show 
how the quantum field theories are related to these classical differential equation theory and 
can be used to obtain some results of them. 

Our conclusion can be summarized as the following points: 

(A) : There are only three asymptotic expansion regions for both the Mathieu and Lame 
eigenvalues, they are in one to one correspondence. 

(B) : The Lame eigenvalue correctly reduces to the Mathieu eigenvalue in the limit n — > 
oo, g — )■ while n'^{n — — )■ h'^/lQ or k — )■ i2h^ keep fixed. 

(C) : The WKB method gives a nontrivial check for the combinatorial approach. This also 
supports the gauge theory/integrable model proposal of [36]. 

(D) : The asymptotic expansion can be consistently extended from the electric region to the 
magnetic and dyonic regions. 

(E) : Eq. (jS]) is more suitable for studying its behavior in the magnetic and dyonic regions, 
while Eq. (|6]) is more suitable for studying its behavior in the electric region. The information 
of dyonic expansion can be obtained from the magnetic expansion by a simple map. 
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8 Appendix: Modular forms, Theta constants 



We collect some basic facts about modular forms and theta constants we have used. Our 
main references are [HI [531 EU ■ 

Theta constant 9i are theta function 9i{z\ q) ai z = 0: 

91 = 0, 

oo 

1 V — ^ n(n + l) 1 Q p. -jj-, ^ ~ 

92 = 2qs2_^q^^=2q^{l + q + q'^ + q^ + q^^ + q^^ + ---), 



n=l 



93 = 1 + 2^9^^ = 1 + 2g^ + + 2g3 +2g^ + 2gf + ■■ ■ , 

n=l 

oo 

^4 = l + 2^(-l)"g'^ = l-2g5 + 2g2_2gi + 2g^-2g¥ + ... . (108) 



n=l 



satisfying 

9t = 9l + 9l 9[ = 929s9^, (109) 
with the nome q = e^'^" , t defined by 

iK' 



K ' 

K and K' are the complete elliptic integrals of the first kind (A; ^1): 



2 



2 



a/1 — A;^ sin^ y9 2 2 2 



(110) 



Jo a/1 — A;^sin^y9 2 2 2 2 

= ![(l + iA;2 + lA;4 + ^A;6 + ...). (Ill) 
2^ 4 64 256 ' ^ ' 



= (l^\e + —k^ + —k^ + ---)\n^ + Oie). (112) 
^ ^ 4 64 256 ^ A; ^ ^ ^ ^ 

where siny^o = |- The complete elliptic integrals of the second kind is 

TT 

E = ^ ' dyp^l - P sin^ = \f{~ ^, 1; fc^) 

= ![(1 _1a;2_ _...). (II3) 

2^ 4 64 256 ^ ^ ^ 
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The relation between the module k, the complementary module k' , and q are 



1 I 

n=l 1 + 2 

= Aq^ - 16g3 + 56q^ - IGOg? + 404g3 - 944g^ + 



n 

4 



oo 



k 



1 - 9"-^ ^4 ^1 



1 ,4 1 ,6 29 ,o 13 11989 ,,3 10879 , 

q = A;^ H k^ H A;^ H A;^° H k^^ H A;^^ H . (114) 

^ 256 256 8192 4096 4194304 4194304 ^ ^ 

The roots Cj = p{zi) of p'^(2;) = 4:p^{z) — g2p{z) — 5^3 = are related to the theta constants 

as 

ei = ^^3^-^^2' = ^ + 16g + 16g' + --- , 

62 = -^(^3^ + ^^) = -^-8g^-8?-32gi-8?2 + ---, 

63 = -^^3^ + ^^|--^ + 8g^-8g + 32gi-8g2 + ... , (115) 
satisfying ei + 62 + 63 = 0. With this we have 

ei — 62 C73 

When we derive expansion of A from we have used 

1 i_l,= _l,^_>_^,»_JlL*.o_JHL,-_... („7) 



61-62 2 32 64 8192 16384 262144 

The second Eisenstein series is represented by 



E^ir) = l-24y ^^ = l-24y , 



n=l n ^ 



= 1 - 24g - 72g2 - 96?^ - 168g^ - 144g^ . (118) 

The equality of the two summation would be obvious from the following, 



00 „ 00 00 00 00 

n=l n=l m=0 n=l m=l 

00 „ ^ 00 00 



E(T^ = E((T^^) = E(E«'")(Eo 

n=l ^ ^ ' n=l ^ ^ ' ^ ' n=l r=l s=0 

00 00 00 00 

n=l r=l s=0 n=l t>l 

the last step uses the fact for a fixed t = r + s, there are exactly t pairs of (r > 1, s > 0) 
satisfy it. E2 is related to the Dedekind Eta function by 

Ai„^(^) = gi?,(^). (120) 
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9 Appendix: Jacobi elliptic functions 

The Jacobi elliptic functions satisfy 



sn^x + cn^x = 1, 



their derivative is 



dn^x+fc^sn^x = 1, 

/c^cn^x+A;^ = dn^x, 
cn^x+fc'^sn^x = dn^x. (121) 



-^snx = cnxdnx, 

CLK 

-— cnx = — snxdnx, 
ax 

-— dnx = — /c^snxcnx. (122) 

Especially, in the WKB analysis we have used 

^-(snxcnxdnx) = 1 — 2(1 + k'^)sD?x + Sfc^sn'^x. (123) 
dx 

The limit of = and k = 1: 

snx|fc=o = sin X, cnx|fc=o = cos x, dnx|fc=o = 1- 



I 1 I 1 

snxL=i = tanh X, cnxL=i = — - — , dnxL=o = cnxL=i = — - — . (124) 

smh X smh x 

10 Appendix: Compare to Langmann's expansion 

E. Langmann developed an algorithm to derive the eigenvalue and eigenfunction of elliptic 
Calogero-Moser(-Sutherland) model for general particle number 115] . In his paper [15] a 
series expansion for the 2-body eigenvalue when g << 1 is given (formulae 26a-d). Although 
his method is different method form combinatorial and WKB methods we present here, we 
will show that his expansion is equivalent to our expansion ( l57|) . Langmann considered the 
2-particle Shrodinger equation 

[-(^ + + ^^(^ " " = (125) 
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with potential V{x) = —p{x) + cq, 

1 °° 1 1 

m=l ^ ^ I 

= ^-2g-6g2-8g3-lV + --- (126) 

After changing variables, the equation becomes the Lame equation describing the relative 
motion. The eigenvalue £ is explicitly given as an expansion in g: £ = £q + £iq + £2(1^ + 
£sq^ + ■ ■ ■ • Especially notice that, in the notation of paper [15], the relation 



2£o = P^ + (m + ns)' (127) 

The term [rii + ^2)^ is the kinetic energy of center-of-mass motion. 
The relation between the eigenvalues £ and B is 

B = -2£ + {ni + + 4:\{\ - l)co. (128) 

if we identify notations a.s P = u, X = n. Then using the relation (11281) to Langmann's 
expansion of we can write our expansion (!57|) in a slightly different form 

3 z/2-1 
8n\n - l)2g2 



(i/2 - 1)3(1.2 _ 4) 

+6(z/2 - 1)2(^,2 _ 2)] 
32n^{n - 1)V 



X [n\n - iy{5u' + 7) - 12n{n - l)(z/^ - 1)^ 



X [2n\n - l)^(9z/^ + 58u^ + 29) 



(i/2 - 1)5(^,2 _ 4) (^2 _ 9) 

-Sn^n - If {u^ - iy{7u^ + 17) 
+An\n - l)2(z/2 - l)2(15z/^ - 371^^ _ 2) 
-24n(n-l)(z/2-l)^(z/2_4) 

+3(z/2-l)4(z/2-3)(z/2-4)] + --- (129) 

Compare fl57p and (11291) . the differences between them is that we have collected some terms 
coming from instanton contribution in f l57|) and they, together with —^n{n — 1)(1 — 2E2{q)), 
equal to |n(n — l)E2{q). Although we have only checked this fact up to order q^, it can be 
checked further by calculating higher order instanton effects in gauge theory. 
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